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Abstract. The construction of effective Hamiltonians arising from Loop Quantum Gravity and 
incorporating Planck scale corrections to the dynamics of photons and spin 1/2 particles is summa- 
rized. The imposition of strict bounds upon some parameters of the model using already existing 
experimental data is also reviewed. 
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m 1. INTRODUCTION 

O ' 

The possibility of bringing quantum gravity induced effects to the observational realm 
has sparkled a lot of attention recently. It is expected that the observation of high energy 
cosmological particles, such as photons OJ] and neutrinos J2] arising from gamma ray 
bursts, will provide the appropriate arena to test such predictions. Also, very precise 
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experiments already performed in atomic and nuclear physics to search for minute 
Lorentz covariance violations have been used to place strict bounds upon such effects 

One of the leading theories providing a consistent description of quantum gravity is 
Loop Quantum Gravity (LQG) 0]. This theory predicts the quantization of space in 
units of £p, with lp being the Planck length [5]. An intuitive way of thinking about this 
is to imagine space being described by discrete cells at very small distances d ~ £p, 
with the standard continuous description being recovered for large distances d > ip. 
From the point of view of a particle immersed in such a space, this granular structure 
will act as an effective media modifying the particle propagation properties with respect 
to those usually assumed in the standard vacuum. Such granularity will induce also 
minute violations of Lorentz covariance, which have been the subject of very precise 
experimental investigations Ja], as well as theoretical descriptions leading to a standard 
model extension which can account for the diversity of observations that have been made 
0]. Modifications arising from LQG constitute a specific realization of such general 
scheme, providing a physical interpretation of the parameters involved. 
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To obtain such modifications starting from the full LQG requires a semiclassical ap- 
proximation where the particles (photons and spin 1/2 particles, for example) are treated 
as classical fields, while an appropriate integration is performed upon the gravitational 
degrees of freedom. In this sense, we are interested in the regime where the matter fields 
are slowly varying while the gravitational variables are rapidly varying. The full Hamil- 
tonian is known in LQG, being a well defined regularized operator acting upon cylin- 
drical functions. These are functions of generalized connections defined upon graphs 
r, characterized by a set of vertices {v} = {vi, V2, • ■ ■ , } and edges {e} = {e 1 , e" ', e'" ', . . .} 
joining those vertices. What is missing is the strict construction of the semiclassical state 
describing the matter field of interest, together with the corresponding large scale con- 
tinuous space-time metric (flat space in our case). With these two ingredients one would 
define and calculate the semiclassical effective Hamiltonian as the expectation value of 
the full LQG Hamiltonian in the corresponding semiclassical state. A rigorous formula- 
tion of this problem has turned out to be complicated and is presently in the process of 
development [8]. Here we take an heuristical point of view, starting from the exact op- 
erator version of LQG and defining its action upon the semiclassical state through some 
plausible requirements. 

Central to our approach is Thiemann's regularization of the LQG Hamiltonian [9]. 
This is based upon a triangulation of space, adapted to the corresponding graphs which 
define a given state. The regularization is provided by the volume operator, with discrete 
eigenvalues arising only from the vertices of the graph. 

The paper is organized as follows 2 : section 2 contains a very compact summary 
of Thiemann's regularization, exemplified in the context of the magnetic sector of 
QED, together with the heuristical scheme employed in our estimations. Section 3 
summarizes the results for the case of photons and spin 1/2 particles and contains a 
brief discussion regarding the choice of some relevant parameters in the model. Section 
4 contains the description of the phenomena from the point of view of the laboratory 
frame attached to earth, moving at a speed v/c ~ 10~ 3 with respect to the Cosmic 
Microwave Background frame (CMB), yielding modifications which can be tested with 
already existing experimental data. 



2. THE CALCULATION 



2.1. The regularized Hamiltonian operator 



The curved space magnetic contribution to the QED Hamiltonian is 
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Limitations of space prevent us to provide a more complete list of references. We apologize to the 
corresponding authors. 



where B a = £ a bcEj, c , F_ bc = d^A^. — d c A b , in standard notation. The underline identifies 
the electromagnetic variables. Here q a b is the three-metric, q = det(q a b) with a,b,c,... 
being spatial indices and Q is the electromagnetic coupling constant. Thiemann's regu- 
larized expression for the above contribution is O 

H B = Y G(v) V ^ KL e MNP x 

zt pU vev{y) v(A)=v(A')=v 



where 
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with being the volume operator. The above equations are written in terms of the 
holonomies hy of the corresponding connections: A l a (gravitational) and A a (electromag- 
netic) along the curves y : .57(A), «/y(A), to be defined below. The triangulation involved 
in © is adapted to the graph T corresponding to the state acted upon, in such a way 
that at each vertex v of T and triplet of edges e,e',e" joining the vertex, a tetrahedron 
is defined with basepoint at the vertex v(A) = v and segments 57(A), / = 1,2,3, in the 
directions of e, e', e" respectively. The arcs connecting the end points of .57(A) and 57(A) 
are denoted ajj(A), so that a loop Ci/j : = si o a\j o sj l can be formed. 



2.2. The semiclassical approximation 

We think of the semiclassical configuration describing the particular matter field (E,B 
in this case) plus flat-space at large distances, as given by an ensemble of graphs T, each 
occurring with probability P(T). To each of such graphs we associate a wave function 
\T,Jf,E,B) = \T,S) which is peaked with respect to the classical electromagnetic 
field configuration together with a flat gravitational metric and a zero value for the 
gravitational connection. In other words, the contribution for each operator inside the 
expectation value are estimated as iflolfTlll 



(r,j?,gB\...q ab ...\r,j?,gi) = dab+ol-^ 

(r,J?,E,B\...A ia ...\r,J?,E,B) = ° + ^(^) > W 

while the expectation values including the electric and magnetic operators are estimated 
through their corresponding classical values E and B. The parameter T > is a real 
number. Not surprisingly, the semiclassical state specifies both the classical coordinate 
and the classical momentum for each pair of canonical variables. The scale J2? >> £p 
of the wave function is such that the continuous flat metric approximation is appropriate 
for distances much larger that jSf , while the granular structure of spacetime becomes rel- 
evant when probing distances smaller that Jzf . Such scale will have a natural realization 
according to each particular physical situation. 



We summarize now the method of calculation [T2I1 . For each graph T 

the effective Hamiltonian is defined as Hr = (r,S|//r|r,S). For a given vertex, inside 
the expectation value, one expands each operator in powers of the segments si (A) plus 
derivatives of the matter fields operators. Schematically, in the case of Q this produces 

H r= I I (r,5|£ W9l (v)...^^..F M (v)f ai ...^---(v,,(A))|r,S). (5) 

vev(r) v(A)=v 

where f contains gravitational operators together with contributions depending on the 
segments of the adapted triangulation in the particular graph. Next, space is considered 
to be divided into boxes, each centered at a fixed point x and with volume Jzf 3 ~ cP x. 
The choice of boxes is the same for all the graphs considered. Each box contains a 
large number of vertices of the semiclassical state (Jz? » £p), but it is considered as 
infinitesimal in the scale where the space can be regarded as continuous. The sum over 
the vertices in is subsequently split as the sum over the vertices in each box, plus 
the sum over boxes. Also, one assumes that the electromagnetic operators are slowly 
varying within a box (J?f « A, with A been the particle wavelength), in such a way that 
for all the vertices inside a given box one can write {T,S\ . . -F_ab( v ) - • • |T, 5) = flF^ix). 
Here is the classical electromagnetic field at the center of the box and \x is a 
dimensionless constant which is determined in such a way that the standard classical 
result in the zeroth order approximation is recovered. Applying the procedure just 
described to © leads to 

H r = IF M (4..(^...Z M P)) 14 1 M' 1+l x 

Box vGBox v(A)=v 

x(r,5|lf ai ...^^-(v,.(A))|r,5), (6) 
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where n + 1 is the total number of factors F pq (x) . The expectation value of the gravita- 
tional contribution is supposed to be a rapidly varying function inside each box. Finally, 
the effective Hamiltonian is defined as an average over the graphs T, i.e. over the adapted 
triangulations : H s = J2 r P(r)Hp. This effectively amounts to average the expectation 
values remaining in each box of the sum ©. We call this average T ai .„ pqpiqi "(x) and 
estimate it by demanding T to be constructed from the flat space tensors 8 a t, and e £ ^ c . In 
this way we are imposing isotropy and rotational invariance on our final Hamiltonian, 
which consequently describes the modified dynamics in a specific reference frame which 
we take to be the CMB frame. Also, the scalings given in © together with the additional 
assumptions: (r,5|...V'...|r,5) — > £p, sf — > tp are used in this estimation. After re- 
placing the summation over boxes by the integral over space, the resulting Hamiltonian 
has the final form 

H S = / d\_F piqi (x)... (>... F M (x)) T a ,.P^-(x). (7) 

Since the approach presented here has made use only of the main features that semi- 
classical states should have, all dimensionless coefficients in the expectation values that 
contribute to T ai _ pqPiqi -(x) in Q remain undetermined. They are subsequently de- 
noted by 0's and ?c's. 



3. THE RESULTS 



Here we summarize the corresponding effective Hamiltonians and modified dispersion 
relations for the cases of photons and two-component spin 1/2 particles. 



3.1. Photons 

A detailed discussion can be found in Refs. | HI 12]. The effective Hamiltonian is 
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The ± signs correspond to the different polarizations of the photon. From the above we 
obtain the speed of the photon ( v±(£,Jzf) = d(d±(k,J>?)/dk ) 



v± = 1 ±40 8 {Up) - 69 3 (k£ P ) 2 + 7 (kip) 
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(10) 



The last expression gives v expanded to leading order in l P , with the estimation Jzf = 
l/k. To first order in (kip) we recover the helicity dependent correction found already in 
the seminal work of Gambini and Pullin lfl3ll . As far as the T dependent terms we have 
either a quadratic (T = 0) or a quartic (T = 1) correction. The only possibility to have a 
first order helicity independent correction amounts to set T = —1/2 which corresponds 
to that of Ellis et. al. [14]. However, we do not have an interpretation for such a value of 
T. 



3.2. Two-component spin 1/2 particles 

The details can be found in Refs. 10. 0]. The effective Hamiltonian is 
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where 
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The corresponding dispersion relation is 
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where the new coefficients are linear combinations of some k's. The velocity 
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within the same approximation. Alternative results based on a string theory inspired 
approach can be found in Ref. B15I1 . 



3.3. The parameters Jz? and T 

In order to produce numerical estimations of some of the effects arising from the 
modifications to the dynamics previously obtained, we must further fix the value of the 
scales Jzf and T. Recall that J?f is a scale indicating the onset distance from where the 
non perturbative states of the spin-network can be approximated by the classical flat 
metric. The propagating particle (photon or neutrino) is characterized by energies which 
probe to distances of order A. In order to preserve the description in terms of a classical 
continuous equation it is necessary that Jzf < X . Two distinguished cases arise: (i) the 
mobile scale, where we take the marginal choice = A and (ii) the universal scale, 
which has been considered in Ref.[ 16] in the context of the GZK anomaly. The study of 
the different reactions involved produces a preferred bound on Jzf : 4.6 x lO^GeV -1 > 
Jzf > 8.3 x 10 _9 GeV _1 . A recent study of the gravitational Cerenkov effect together 
with neutrino oscillations ifT^l produces a universal scale estimation which is consistent 
with the former . Bounds for T have been estimated in Ref. Ill] based on the observation 



that atmospheric neutrino oscillations at average energies of the order 10~ 2 — 10 2 GeV 
are dominated by the corresponding mass differences via the oscillation length L m . This 
means that additional contributions to the oscillation length, in particular the quantum 
gravity correction Lqg, should satisfy Lqq > L m . This is used to set a lower bound upon 
T. Within the proposed two different ways of estimating the scale jSf of the process we 
obtain: (i) T > 0.15 when Jz? is considered as a mobile scale and (ii) 1.2 < T when the 
scale Jzf takes the universal value Jzf « 10~ 8 GeV 1 . 



The previously found Hamiltonians were obtained under the assumption of flat space 
isotropy so that they account for the dynamics in a preferred reference frame. We 
have identified it as the frame in which the Cosmic Microwave Background looks 
isotropic. Our velocity w with respect to that frame has already been determined to be 
w/c ~ 1.23 x 10~ 3 by COBE. Thus, in the earth reference frame one expects the appear- 
ance of signals indicating minute violations of space isotropy encoded in w-dependent 
terms appearing in the transformed Hamiltonian or Lagrangian B3J]. On the other hand, 
many high precision experimental test of rotational symmetry, using atomic and nuclear 
system, have been already reported in the literature. Amazingly such precision is al- 
ready enough to set very stringent bounds on some of the parameters arising from the 
quantum gravity corrections. In Ref. |0] we have considered the case of non-relativistic 
Dirac particles obtaining corrections which involve the coupling of the spin to the CMB 
velocity together with a quadrupolar anisotropy of the inertial mass. The calculation was 
made with the choices T = and = 1 /M, where M is the rest mass of the fermion. 
Keeping only terms linear in ip, the equation of motion arising from the two-component 
Hamiltonian (fTTT) can be readily extended to the Dirac case as 



where we have used the representation in which 75 is diagonal, the spin operator is L = 
(i/2)en m y l Y n , K = ®$m 2 lp and a = 1 + &3in£p. The normalization has been chosen 
so that in the limit (m£p) — > we recover the standard massive Dirac equation. The term 
m ( 1 + 03 mlp) can be interpreted as a renormalization of the mass whose physical value 
is taken to be M = m ( 1 + ©3 mtp) . After this modification the corresponding effective 
Lagrangian is 



which describes the time evolution as seen in the CMB frame. In order to obtain the 
dynamics in the laboratory frame we implement an observer Lorentz transformation. To 
this end we rewrite (fToT) in a covariant looking form, by introducing explicitly the CMB 



4. OBSERVATIONAL BOUNDS USING EXISTING DATA 
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frame's four velocity = y(l, w/c). In the metric with signature —2 the result is 
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From Eq. (26) of we obtain the non-relativistic limit of the Hamiltonian correspond- 
ing to (fT71) . up to first order in i P and up to order (w)/c 2 , which is 
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where s = a/2. Here we have not written the terms linear in the momentum since they 
average to zero. In (fT8l) g is the standard gyromagnetic factor, and Qp is the momentum 
quadrupole tensor with components Qpjj = pipj — l/3p 2 8ij. The terms in the second 
square bracket represent a coupling of the spin to the velocity with respect to the "rest" 
(privileged) frame. The first one has been measured with high accuracy in references |@] 
where an upper bound for the coefficient has been found. The second term is a small 
anisotropy contribution and can be neglected. Thus we find the correction 
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The last term of (fTSl) . which represents an anisotropy of the inertial mass, has been 
bounded in Hughes-Drever like experiments. With the approximation Qp = — 5/3 < 
p 2 > Q/R 2 for the momentum quadrupole moment, with Q being the electric quadrupole 
moment and R the nuclear radius, we obtain 
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for the quadrupole mass perturbation, where 9 is the angle between the quantization axis 
and w. Using < p 2 /2M >~ 40 MeV for the energy of a nucleon in the last shell of a 
typical heavy nucleus, together with the experimental bounds of references Jal we find 



02+^04 
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(21) 



The above bounds on terms that were formerly expected to be of order unity, already 
call into question the scenarios inspired on the various approaches to quantum gravity, 
suggesting the existence of Lorentz violating Lagrangian corrections which are linear in 
Planck's length. To this respect it is interesting to notice that a very reasonable fit to the 
gamma ray spectrum beyond de GZK cutoff has been recently made by using dispersion 
relations of higher order than linear in £p [|19il . Observational bounds upon parameters 
of related theories are obtained in the works of Ref. Il20ll 
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